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AN  INFORMAL  INTRODUCTION 

The  positioning  problaa  arises  when  It  is  necessary 
to  locate  a sat  of  geographically-distributed 
objects  using  aeaaureaents  of  the  distances  between 
soma  oo j eat  pairs.  In  a Packet  Radio  Network,  for 
instance,  any  two  network  members  that  can  talk  to 
aeon  other  may  use  a staple  time-stamping  aechanlso 
to  measure  the  distanae  between  them;  a distance 
measurement  protocol  say  then  be  developed.  The 
problem  is  whether  and  how  the  distanae 
measurements  can  be  used  to  determine  the 
geographical  location  with  respect  to  a given 
system  of  coordinates. 

A knowledge  of  the  precise  loaation  of  eaah  network 
node  Is  crucial  to  the  operation  of  Distributed 
Sensors  Networks.  The  data  oolleoted  and 

interpreted  by  different  sensors  may  be  correlated 
and  integrated  only  If  we  know  their  preoise 
looatlon.  A posit ion- locating  syatea  may  be 

Invaluable  to  the  operation  of  a fleat  of  vehicles, 
each  equipped  with  a Paoket  Radio  Unit.  For 
maple,  totiitoring  the  looatlon  of  a fleet  of 
security  vehicles,  aircraft,  a tank  division,  or  a 
flock  of  missiles  could  all  be  assisted  by  a 
position-locating  syatea.  Clearly  a positioning 

system  would  be  in  Important  service  to  Packet 
Radio  Network  users.  v , 

A few  proolera  must  be  solved  before  a good 
positioning  system  say  be  developed: 

l.  efficient  algorithms  to  determine  the 
location  of  objects  by  using  distance 
measurements  would  be  developed. 

I.  Conditions  under  which  a solution  exists  or 
does  not  exist  snould  be  identified. 

I.  Conditions  under  which  there  exists  e 

unique  solution  should  be  established. 

4.  Conditions  under  which  there  exists  s 

finite  nuabsr  sf  solutions  snould  be 
Identified.  It  should  also  be  understood 
how  to  transfora  ona  solution  Into  another. 

3.  Conditions  under  which  the  solution  is 

insensitive  to  saall  measurement  errors 
snould  be  established. 

6.  Tight  bounds  upon  the  acaurtcy  of  the 
solution  should  be  determined. 

?.  Ill-conditioned  problems  should  be 
Identified. 


However,  while  the  formulation  of  tbe  problems  Is 
simple,  the  mathematical  and  algorithmic 
intrioaoles  of  deriving  solutions  are  perplexing. 

To  develop  some  insignt  into  the  problems,  let  us 
consider  a few  simple  examples.  Tbe  simplest 
positioning  problem  of  interest  Is  to  looate  three 
points  using  distanae  measurements.  By  means  of 
simple  trigonometry,  this  may  usually  be  done 
easily.  However,  let  us  consider  a degenerate 
triangle  (Figure  1).  Because  the  system  is  very 
sensitive  to  errors,  a small  error  in  the 
meesureaent  may  produoe  a large  error  in  the 
computed  position.  Some  of  the  questions  to  be 
addressed  are  as  follows:  Why  Is  tha  degenerate 
triangle  sensitive  to  errors?  How  can  we  determine 
whether  or  not  other  systems  are  sensitive? 


Figure  1 . A Degenerate  Triangular  Syatea 

Positioning  systems  may  be  constructed  by  a siaple 
prooadure  of  pasting  triangles  together,  and  auah 
systems  may  be  positioned  by  solving  the  triangles 
from  which  they  are  constructed.  For  instance, 
consider  the  system  of  points  depicted  In  Figure  a. 
It  Is  possible  to  Looate  the  points  In  ths  order 
numbered.  However,  the  seme  systea  admits  a few 
solutions  (the  number  of  wnich  grows  exponentially 
with  the  >mber  of  nodes) . If  we  had  some  further 
Information  scout  the  positions  of  the  objects,  now 
could  we  use  It  to  Identify  the  true  solution?  For 
instance,  If  one  node  is  Unown  to  be  a vehicle 
moving  on  a certain  road,  sany  of  the  feasible 
solutions  that  satisfy  the  Slstanoe  constraints  can 
be  eliminated,  because  they  assign  the  vehicle  to  a 
position  not  on  tha  road.  How  should  that 
elimination  be  effected? 


Figure  e.  A Triangulated  Positioning  Syatea 


U? 


In  particular,  If  a now  measurement  ware  given 
(Figure  3),  now  can  wa  find  the  right  solution 
among  ail  possible  ones?  (A  few  solutions  any  match 
the  measurements  to  within  a given  error.)  Me  have 
shown  that  this  last  decision  problem  belongs  to 
the  olass  of  difficult  combinatorial  decision 
problems— -the  so-called  NP  complete  problems.  The 
last  statement  also  proves  that  the  existence 
prooiea  (i.a..  Does  a solution  exist?)  is  NP 
complete. 


Figure  3.  A triangulated  Positioning  System  with 
an  additional  Measurement 


contains  a subsystem  that  may  be  solved 
independently.  Namely,  any  incremental 
construction  algorithm  must  decide  whether  or  not  a 
given  positioning  system  is  primitive  (i.e., 
constructible  but  having  no  construe  tide 
subsystems).  the  last  problem  seems  to  be  a 
difficult  combinatorial  problem  which  we  suspect  to 
be  NP  complete  (though  we  do  not  yet  wnou  how  to 
prove  this). 

If  the  last  conjecture  is  true,  then  the  problem  of 
constructing  solutions  to  t.he  positioning  problem, 
using  an  incremental  algorithm,  is  NP  complete. 
This  unfortunate  result  does  not  Imply  that  "brute 
foroe*  (i.e.,  iteraf  've  algorithms)  should  be 
preferred.  It  is  reasonable  to  believe  tnat  most 
actual  positioning  problems  may  be  better  solved  by 
means  of  an  Intelligent  incremental  algorithm.  The 
precise  meaning  of  "most"  is  yet  to  be  defined. 

Numerous  other  challenging  and  interesting  related 
problems  exist.  While  we  will  not  mawe  a 
comprehensive  presentation,  we  will  examine  some  cf 
the  problems  formally,  expose  the  difficulties,  and 
present  some  partial  solutions  we  have  developed. 
This  report  is  an  extended  summary  of  our  present 
state  of  Knowledge.  A more  detailed  report  la  now 
being  prepared. 

1.  THE  PROBLEMS 


Not  ill  positioning  ajdems  may  be  solved  with  the 
tid  cf  triangles.  In  r*ot,  for  a system  possessing 
* sufficiently  large  number  of  nodes,  it  is  always 
possible  that  the  position  of  the  nodes  can  be 
located  only  by  solving  for  the  location  of  all 
points  simultaneously.  Suoh  unfortunate  systems 
require  an  enormous  amount  of  computation.  For  a 
simple  example,  consider  the  hexagon  af  Figure  >». 
It  is  Impossible  to  solve  the  location  of  Its  nodes 
using  an  incremental  algorithm;  all  must  be  solved 
at  once. 


Figure  ».  A Positioning  System  which  say  wily 
se  solved  simultaneously 


fortunately  enough,  such  primitive  systems  (whose 
parts  aannot  be  positioned  unless  the  whole  system 
is)  seem  to  be  rare.  Many  positioning  systems 
sould  se  solved  using  an  incremental  process  (wnicn 
simplifies  the  solution  algorithm  and  increases  its 
speed  ana  accuraoy).  However,  an  algorithm  that 
would  construct  the  location  of  a given  point 
system  by  ionstruotlng  subsystems  first  would  be 
able  to  identify  constructiDle  parts.  In 
partloular.  suen  an  algorithm  should  ee  able  to 
decide  whether  or  not  a given  positioning  system 


The  positioning  problem  can  be  described  as 
follows; 

i.  p =S  (Pr  P,,  . . . Pj,],  e sat  or  points  in  • 
the  plane. 


2. 


A set  of  distance  measurements  between  some 
pairs  of  points.  Each  measurement  datum 
consists  of  the  ideatity  of  the  peir  ?,  and 


?j,  the  measured  distance  d ^ and 
estimate  of  the  measurement  error  « 


an 


U* 


3.  Position  coordinates  for  at  least  three 
points,  say  P,,  Pj,  Pj,  to  be  called  the 

mi  ir.uaa.la. 


Me  snail  call  the  set  af  points  ?,  together  with 
the  distance  measurements  Id,,)  and  the  base 
triangle,  a acini  system. 

1.2  PSCELcHS 

a aaalUU  aaaiuca  af  system  is  » set  ?r 

coordinates  that  satisfies  the  dlstanae  constraints 
and  assigns  to  the  base  triangle  its  actual 
o me rdi nates.  The  set  of  ail  feasible  positions 

will  be  called  the  splupjan  j£L  of  the  positioning 
problem.  The  positioning  problem  consists  of 
chareo ter icing  the  solution  set  namely , 

1.  Is  the  solution  set  ebpty?  ( ax intones? 

2.  hoes  the  solution  set  contain  a eor.ti.nuua 
of  solutions,  or  is  it  * dlsorete  set)  la 
U s finite  set?  UcnsrsLlsea.  unitucr.sss) 
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3.  If  the  solution  sat  is  finite,  what  are  all 
the  feasible  solutions?  (position 

aaaatrusUaa) 

4.  What  is  the  error  in  the  position  due  to 
propagation  of  measurement  errors?  (error 
analysis) 

5.  What  is  the  sensitivity  of  the  solution  set 
to  measurement  errors?  ( sensitivity 
analysis) 

In  what  follows  we  examine  some  partial  answers  to 
some  of  the  difficult  problems  posed  above. 

Note  that  in  the  sequel  we  restrict  ourselves  to 
the  problem  of  oositioalng  points  relative  to  each 
other,  i.a.,  with  respect  to  any  coordinate  system 
of  our  onoice.  This  leaves  us  three  degrees  of 
freedom  (two  for  translation  and  one  for  rotation) 
and  the  orientation  for  our  choice  of  the 
coordinate  system.  With  the  aid  of  the  third  item 
on  the  input  list  (section  1.1)  it  is  possible  to 
position  tne  point  system  absolutely  once  It  has 
been  positioned  relative  to  an  arbitrary  coordinate 
system. 

2.  GEOMETRIC  RESULTS 

We  associate  with  the  point  system  a graph  whose 
vertiaea  represent  points  and  whose  edges  represent 
distance  measurements.  We  call  this  graph  the 
aaasuraawnta  graph,  and  say  that  a given  property  of 
a point  system  is  ssBDlBttgrial  when  it  asn  be 
expressed  In  taros  of  the  measurements  graph  only. 
A structure  is  a graph  together  with  a sapping  of 
edges  into  positive  real  numbers  which  u*  call 
lengths.  A positioning  system  say  oe  considered  as 
a pin- Jointed  bar  structure,  t.a.,  a truss. 
Problems  of  uniqueness  of  (l.e.,  structure  of  the 
solution  set)  for  the  positioning  problem  translate 
Into  problests  of  rigidity  of  the  respective  truss. 
Problems  of  a solution's  sensitivity  to  errors 
trsnslats  into  problems  of  infinitesimal  rigidity 
of  cite  respective  truss  (l.e.,  admissibility  of 
Infinitesimal  flexing  of  the  truss).  Problems  of 
constructing  a solution  to  the  positioning  problem 
correspond  to  aonstruotlen  of  the  truss.  Therefore 
we  snail  use  methods  and  terminology  that  pertain 
to  both  structures  and  trusses. 

2.1  RIGIDITY 

The  results  In  this  area  fall  Into  three  classes: 

1.  1 number  of  different  oharadterlsatlons  of 
infinitesimal  rigidity  (error  sensitivity). 

2.  An  algorithm  to  determine  whether  or  not  e 
given  solution  of  the  positioning  problem 
is  sensitive  to  measurement  errors. 

3.  A tcsbinatorlel  cnarioterltatlan  of  plant 
rigidity  m teras  of  a property  of  the 
underlying  aeasureeants  graph. 

While  the  pro-nea  of  structural  rigidity  has 
attracted  mathematicians,  er slr.ee rs,  and  arcnitects 


for  several  centuries  A the  solution  to  most 
fundamental  questions  of  rigidity  are  far  from 
tcnowu.  Reoently  interest  in  this  age-old  problem 
has  been  renewed  (WHITELT  77,78],  and  some 

significant  contributions  produced. 

2.1.1  Characterization  of  rigidity 

The  first  problem,  l.e.,  that  of  characterising 
rigidity,  has  a few  solutions,  all  of  which  employ 
local  infinitesimal  characterizations . Loosely 
speaking,  a structure  is  rigid  if  it  does  not  admit 
relative  motions  of  its  parts,  that  is,  if  the  only 
motions  which  it  admits  are  trivial  (l.e.,. 

translations  and  rotations).  Therefore  the  study 
of  rigidity  is  a study  of  possible  motions.  A 
rigid  structure  corresponds  to  a positioning 
problem  with  a discrete  solution  set.  An 
infinitesimally  rigid  structure  corresponds  to  an 
error-insensitive  solution. 

There  are  two  approaches  to  motions  of  structures: 
It  is  possible  to  consider  the  velooity  veotors  of 
the  nodes  or  the  relative  angular  motion  of  edges 
attached  to  a common  node.  Accordingly,  it  is 
possible  to  develop  two  notions  of  infinitesimal 
rigidity.  Another  possible  approaan  Is  to  consider 
the  stresses  in  the  structure  resulting  from 
applying  external  forces.  Rigidity  may  be  defined 
es  the  ability  of  the  structure  to  resolve  forces. 
It  is  possible  to  show  that  both  the  aoove 

approaches  are  equivalent  (GLUCS  75,  UHITELY  77, 
73]. 

OOSG  nmbleas: 

1.  Characterization  of  rigid  structures  which 
are  not  Infinitesimally  rigid. 

An  infinitesimally  rigid  structure  1 Is 
rigid,  out  not  vice  versa.  It  is  possible 
for  a structure  to  be  rigid  (l.e,,  admit  no 
finite  relative  motion*  of  its  parts)  yet 
to  admit  Infinitesimal  perturbations  (l.e., 
be  sensitive  to  errors).  A typical  example 
of  in  error-rigid  rigid  structure  is  tne 
degenerate  triangle  la  Figure  1. 

The  difficulty  in  solving  this  preoles  is 
that  we  do  not  possess  extensive  tools  far 
global  analysis,  while  local  analysis  is 
wall  developed. 

2.  Characterisation  of  rigidity  with  respect 
to  discontinuous  motions  such  as 
reflections. 

Consider  tht  pair  of  pasted  trlsngles  m 
Figure  3(a)  salow.  The  two  trlanglas  say 
be  positioned  with  respect  to  each  other  m 
two  distinct  ways.  The  two  resulting 
structures  are  rigid  ‘do  not  adait  non 
trivial  actions;  but  admit  relative 
reflection  of  parts.  Other  structures 


^ A partial  list  of  resaarohers  interests  In  the 
problem  includes  Pascal,  Suler,  Causny,  Xasueil, 
Cayley . Alexandrov,  and  others. 
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admit  a more  aonplax  fora  of  discrete 
mo  van  aat  of  their  internal  parts.,  e.g., 
Figure  5(b) . 

the  problem  of  characterising  rigidity  with 
respect  to  discrete  motions  Is  difficult, 
for  we  cot  only  have  to  address  a problem 
of  global  analysis,  but  also  face  difficult 
problems  of  combinatorial  topology. 


4 3 


3(b) 


Figures  3(e)  end  (b).  Two  cases  of  two  rigid  struc- 
tures solving  the  sue  positioning  problem. 


J.i.O  error  sensitivity 

Ute  different  definitions  of  infinitesimal  rigidity 
induoe  different  algorithms  to  determine  unether  or 
not  a given  structure  is  lnflnitesiasily  rigid. 
Houever,  most  of  those  algorithms  have  been 
produced  by  and  for  math  emit  ideas  unconcerned  with 
tensputat  lonal  efficiency.  We  nave  developed  a 
novel  algorithm  to  teat  unether  a given  feasible 
solution  of  the  positioning  problem  Is 
infinitesimally  rigid,  l.m.,  insensitive  to  errors. 

The  idee  behind  the  rigidity  testing  algorithm  is 
simple)  try  to  solve  for  an  admissible  assignment 
of  infinitesimal  velocities  that  fle*«a  the 
structure.  It  is  necessary  to  examine  only  the 
effects  of  e velocity  assignments  over  e set  of 
basic  circuits  of  tbe  underlying  measurement  grepB 
in  order  to  reduce  the  problem  to  the  solubles  of  a 
Llneai'  system  of  equations. 


Oran  groLlras; 


1.  SstibUsh  measures  of  error  sensitivity  and 
algorithms  to  confute  sensitivity. 


The  characterisation  of  error  sensitivity 
In  terns  'of  infinitesimal  rigidity  is  too 
crude.  We  would  nice  to  have  an  estimate 
of  iiax  much  error  sensitivity  a given 
structure  possesses. 

2.  Develop  sensitivity  measures  in  terns  of 
the  distance  measurements  data,  not  the 
particular  solution  they  yield. 

In  addition  to  the  difficulty  Involved  in 
developing  apriorl  sensitivity  measures,  we 
face  the  difficulty  that  the  siae 
measurement  data  may  have  a number  of 
solutions,  each  possessing  different 
sensitivity.  Apriorl  it  is  even  possible 
for  a given  system  of  measurement  data  to 
possess  some  rigid  and  some  nonrlgld 
solutions.  Figure  6 depiots  tuo  solutions 
of  the  sane  positioning  problem,  one  rigid 
and  the  other  infinitesimally  flexible. 


t 


Figure  6.  1 positioning  problem  possessing  a rigid 
and  an  lofinltMimaU y-desiblo  loiutioa 


We  snail  now  dmmorltes  the  third  alas*  of  results  in 
acre  detail. 

2.2.  ?SGH  SfiCKSTSS  TO  CCMSISIATOSICj 
2.2.1  Stiff  graph* 

The  so st  powerful  results  of  t»s  xeudy  of  rigidity 
appear  to  be  combinatorial  cnaraotertsattona  of 
rigid  structures.  T»e  ideas  behind  tne  passage 
from  geometry  to  combinatorics  are  founded  :n  scat 
simple  intuitive  experiences ? the  major  idee  vf 
union  is  mat  wae  cosblnauons  of  bars  (edges)  and 
hinge*  (nodes)  ere  rigid  for  almost  any  choice  of 
plane  Labeddtng.  For  instance,  the  fuu  graph  en 
three  nodes  is  rigid  for  ell  plan#  Sabeddlngs;  It 
is  not  inf Initealsally  rigid  uhefi  tbe  three  points 
sre  eo linear . Slallarl*,  structures  uncs# 


intuition  should  net  be  pursued  blindly,  however, 
unsa  It  teats  to  prealems  cf  rigidity,  fer—as  we 
snail  see — many  intuitive  expectations  turn  cut, 
surprisingly,  to  be  false. 


1W 


.■  ■, . r. » . 
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underlying  graph  la  a full  graph  are  rigid  for 
almost  any  piano  imbedding.  On  tha  other  hand, 
soma  graphs  will  produce  flexible  structures  for 
almost  any  plane  imbedding  (e.g.,  a circuit  on  four 
nodes  is  alaoat  always  flexible,  except  whan  one  of 
the  edges  is  assigned  a zero  length  or  when  all 
four  points  are  colinear).  The  problem  is:  la  li. 
possible  la  nharnntertse  granha  almost  ail  Si  whose 
emhwddings  fora  rigid  structures?.  This  question 
has  been  the  major  problem  of  the  theory  of 
rigidity. 

We  shall  define  a critical  combinatorial  property 
of  graphs  which  we  call  (plana)  stiffness.  First 
we  associate  with  a graph  3 £ <V,E>  a nuaber 

f(3)  3 2 1 V I - 1 El -3 , measuring  tha  ovarall  excess 

of  unknowns  over  constraints.  (Each  vertex 
contributes  two  unknown  coordinates,  and  each  edga 
constrains  tha  two  vertices  inoident  unto  it 
through  a single  equation  for  the  distanoe.  Note 
that  we  discount  three  degrees  of  freedoa  to 
account  for  possible  external  motion,  i.a., 
translation  and  rotation.)  The  quantity  f(G) 
measures  the  overall  freedoa  if  Internal  movement 
of  the  graph. 

Tha  quantity  f(3)  may  be  used  to  express  the 
property  of  stiffness,  i.a.,  of  having  a sufficient 
dumber  of  constraints  to  prevent  relative  motions 
of  different  parts  of  a graph,  loosely  speaking,  a 
graph  is  stiff  if  it  is  poseible  to  remove  soae 
redundant  edges  so  that  the  remaining  graph  has  0 
dejreas  of  internal  freed:*,  and  non  of  its 

su'Otfraphb  has  an  exctsa  of  constraints  Cl. a.,  a 
negat iv«  internal  freedom).  Formally,  a graph 
Oe<V,S>  is  st^ff  if  it  has  a spanning  subgraph 
■VaCV.K’V  Cl.*,.  3'  is  gaoarated  by  ramoviag 

excessive  constraints  fro*  3)  suoh  that 

l.  f'iOMtQ  (i.a.,  G*  mu  0 degrees  of  internal 
freedom) 

i.  if  3"  la  a ay  subgraph  of  3'  then 

f(iS*)  £ 0.  Cl. a.,  3*  does  not  possess 
internally  ®v*r-ccnst rai ued  subgraphs) . 


3. i.J  The  rigidity  toewea 

The  meat  important  reault  of  the  geometric  theory 
of  posit leala#  is 

ti&QSEH  CFlaae  ill  gutty!  J: 

k tu&A  » [tMa  ntxft  * pun 
stmsam&JL  maiu  fomrulx*.  alacai 
ill  im  ixtaltim.  at  a mil 
asm 

(Hare  *ai»3St  ill*  is  used  in  the  topological 
sense,  i.«.,  the  set  of  rigid  plane  isbeddlngs  of  * 
'•Iff  graph  is  open  and  ieaw  in  the  space  of 
imbedding*,  which  further  implies  that  for  any 
isrei  probability  aeasure  on  the  space  of  all 
iaosddings,  ooatinuaus  with  respect  to  Lebsague 
seaaure,  the  set  of  «an  rigid  Isbeddlngs  of  a stiff 
graph  is  of  measur « >i.) 


The  theorem  above  is  a signifioant  tool  for 
handling  positioning  problems,  namely,  It  makes  it 
possible  to  infer  properties  of  structure  and 
derive  answers  to  the  problems  of  positioning  by 
examining  tha  measurement  graph  only.  The  results 
that  we  derive  will  he  true  for  almost  all 
assignments  of  distacaes  to  edges  of  the 
measurement  graph,  which  greatly  simplifies  the 
study  of  the  positioning  problem. 

To  further  appreciate  the  power  of  the  result  let 
us  note  in  passing  that  the  theorem  does  not 
generalize  even  to  ’three  dimensional  structures. 
That  is,  it  is  possible  to  have  a sufficient  number 
of  well  distributed  constraints  and  yet  haw  a 
structure  with  a continuum  of  solutions  no  matter 
what  lengths  are  assigned  to  the  edges.  Figure  7 
below  depicts  a typical  system.  The  combinatorial 
characterization  of  rigid  structures  in  spaaas  of 
dimensions  greater  than  two  la  an  open  problem. 


Figure  7.  k counter  example  to  a 3 -dimensional 
rigidity  theorem 


3.3  LIMITATIONS  CP  THE  RIG ID ITT  THEOREM 

The  rigidity  theorem  guarantees  that  a structure 
baaed  on  a stiff  graph  will  aieost  always  be  rigid. 
However,  it  is  possible  to  assign  lengths  tc  edges 
of  a stiff  graph  such  that  tha  resulting  structure 
admits  iofinitealaal  flexing  and  is,  tharefors 
vary  sensitive  to  errors.  In  fact,  the  structure 
becomes  so  error -increasing  seohanlsa.  One  such 
structure  is  the  degenerate  triangle  of  Figure  1 , 
another  is  Pascal's  hexagon  depicted  in  Figure  3. 
This  hexagon  is  stiff  and  thus  rigid  for  almost  ail 
plane  labeddings.  However,  the  hexagon  is 
infinitesimally  flexible  whenever  (and  only  if)  its 
nodes  lie  on  a conic  section.  This  dinar  re  result, 
due  to  [CROFTGN  t873],  has  been  rediscovered 
Independently  by  cany  researchers  (including 

ourselves).  The  preef  follows  frea  a simple 

application  cf  a celebrated  theorem  of  Pascal. 

in  even  worse  esse  is  that  cf  a stiff  graph 
admitting  a continue  notion  in  imbedding 

of  the  stiff  graph  that  is  net  even  gioSaUv 
rigid).  Such  a structure  is  depicted  in  Figure  ?. 
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Figure  9.  4 flexible  structure  whose 
underlying  graph  la  rigid. 
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llthougn  the  underlying  graph  la  stiff  (even 
overconatrained) , the  solution  set  of  the 
respective  positioning  probleei  contains  a continuua 
of  solutions. 

To  produce  suon  an  example  we  started  with  * 
.•wnrigld  structure  (l.e<»  a aeahsnisa),  than 
setlcuiously  added  Sara  that  did  not  conatraia  the 
notion  of  tne  original  seehanlsm.  Specifically,  we 
started  with  a four-bar  aeentnia*,  selected  a point 
on  any  diagonal  and  connected  this  point  to  points 
an  tha  original  se  eft  anise.  3y  careful  selection  of 
the  aonneatlona  to  produce  two  parallel ogress,  the 
action  of  the  original  aechanlse  la  not  perturhed 
oy  the  additional  ‘constraints.*  Iddltlanal  sera 
are  uaed  only  to  hold  each  of  the  original  four 
Cara  together. 

The  rigidity  theoreo  guaranteed  that  such 
'unpleasant  poattloaing  aystaes  are  extremely  rare, 
let  we  should  tear  la  alnd  that  they  say  exist 
(highly  syasatrlo  structural  are  very  likely 
exceptional . 

}.  CCMBXSAtORICS  3F  PGSlTSOMXKil 

Having  aeen  the  algclflcacca  and  limitations  of 
atlffneaa,  we  will  now  atudy  stiff  graphs  In  order 
to  develop  aethcds  for  recognising  stiffness  and 
for  tearing  structures  into  stiff  Suiparta,  later 
cesantlog  those  parts  together.  (Tearing  Is  tha 


natural  instrument  for  a dlvide-aad-conquer 
approach  to  the  construction  problem.) 

3.1  ciuiUdEaixmoH  of  stiff  oa*p«s 

The  question  that  we  would  like  to  address  in  this 
Motion  Is:  Iflum  1 graur  stiff?.  We  wish  to 

dertva  nsoessary  and  suffloiir.t  aonditlans  for 
stiffness  in  terms  of  well  known  grapn  properties. 
The  cost  natural  of  classical  grapn  properties 
relating  to  stiffness  Is  connectivity . 

3.1.1  3om  slmle  cearsoteri  cations 

We  have  Seen  ahle  to  derive  e list  of  useful 
properties  of  stiff  graphs: 

1.  4 stiff  graph  is  a clock  tl.s.,  has  no  cut 
vertex) . 

2.  4 out-set  of  a stiff  gfsph,  separating  it 
Into  two  nontrivial  coapononts,  must 
contain  at  least  3 edges. 

3.  4 3-out-eet  of  a stiff  graph,  separating  it 
into  two  nontrivial  icaponents,  oust  also 
separate  it  Into  two  stiff  graphs. 

c.  if  3«<V,S>  is  a stiff  grapn  and  V*  s V ■ 
vertex  out-set  separating  3 into  tcapcnents 
1,  • <Vl,£l>  and  if  the  subgraph  of  3 
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spanned  by  V'  la  stiff,  than  ao  am  the 
subgraphs  apannad  by  V'  u Vi. 

5.  If  v la  a vartex  of  degree  two  in  a graph 
0,  than  3 la  atlff  iff  the  subgraph  of  G 
formed  by  removing  v is  stiff  too. 

The  above  results  are  a sample  from  a larger  class 
of  results,  all  of  which  serve  to  establish  tools 
for  a "divide  and  conquer"  approach  to  the 

stiffness  problem.  For  Instance,  we  would  like  to 
determine  whether  a given  graph  can  be  torn  into 
stiff  aubparta  wtiiob  may  later  be  used  to 

synthesize  the  original  graph.  Figure  10  depicts  a 
typical  configuration,  corresponding  to  result  (3) 
in  the  above  list.  Figure  11  depicts  another 
possible  configuration  which  admits  tearing,  this 
time  a particular  instance  of  (b). 

3.1.2  Stiffness  and  connectivity 

Loosely  spoaxinj,  stiffness  is  a property  of  graphs 
wnicn  has  to  do  with  the  density  of  edges,  i.e.,  it 
is  a measure  of  haw  well  different  nodes  ere 
ittaohed  to  each  other.  It  is  only  natural  to 
expect  that  such  a property  should  hear  relation  to 
classical  measures  of  edge-density. 

There  ere  three  major  classical  measures  of 
edge-density:  node  degrees,  minimal  edge  cut-set, 

minimal  vertex  cut-set.  Here  we  explore  the 
relations  eaong  these  three  properties  and 

atiffneas. 


1.  Thera  exist  graons  whose  nodes  possess 
arbitrarily  large  degrees  but  which  are  not  stiff. 

To  prove  this  result  ue  describe  a simple  method  to 
construct  counter  examples.  Start  with  a flexible 
graph,  say  a four-bar  mechanism.  Add  nodes  and 
edges  to  increase  the  degrees,  preserving  the 
flexibility.  This  process  is  demonstrated  ia 
Figure  12. 

The  above  process  serves  to  show  that: 

2.  There  exist  graphs  with  aa  arbitrarily  large 
minimal  edge-out-set  but  which  are  not,  stiff. 

Thus,  two  of  the  classical  measures  of  connectivity 
are  not  related  to  stiffness.  Let  us  now  consider 
the  strongest  measure  of  connectivity, 
vertex-conneotivi ty . 

We  have  seen  in  the  previous  seatlon  that  a stiff 
graph  la  at  least  2-aonnacted.  It  Is  easy  to 
construct  2 and  3-conneoted  graphs  wnich  are  not 
stiff.  But  with  ic-conneoted  graphs,  '<  > b,  the 
problem  is  ao  longer  simple.  Indeed  any 
b-connected  graph  must  have  a minimal  degree  which 
is  at  least  4.  thus  the  total  number  of  edges  a la 
at  least  twice  the  number  of  nodes.  Therefore  the 
overall  number  of  internal  degrees  of  freedom 
f(G)a2n-«-3  is  not  greater  than  -3-  Hot  only  Is  a 
b— connected  graph  over-constrained,  but  the 


figure  '1.  Tearing  alcng  t stiff  vertex  cut-set 
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connectivity  impliss  that  the  edges  must  he  wall 
distributed.  Intuitively  one  would  expect  that  a 
4-conneoted  graph  is  stiff. 

Not  so.  It  is  possible  to  construct  4-conneoted 
and  even  5-conaeoted  graphs  which  are  flexible. 
Two  such  examples  are  depicted  in  figures  13  and 


The  process  we  applied  to  derive  these  two 
surprising  graphs  cannot  be  applied  to  produoe  5- 
(or  sore)  conneoted  graphs  whioh  are  flexible.  Vo 
do  not  know  at  all  unether  such  graphs  even  exist. 
The  problem  is  open: 

3.  Is  there  a number  It  such  that  any  k-oonnwcted 
■itraoh  la  stiff? 

It  this  point,  however,  the  problem  is  mainly  of  an 
academic  interest,  for  a condition  which  requires 
such  a high  connectivity  seems  to  be  of  ho 
practical  significance. 

To  summarise,  we  have  seen  that  the  relation 
cetween  stiffness  and  measures  of  connectivity  (if 
there  is  any)  are  not  simpie,  contrary  to  the 
ipnorl  intuition  which  leads  us  to  explore  these 
relations. 


3.2  CONSTRUCTION  PROBLEMS 

i position-locating  algorithm  is  essentially  a 
process  tnat  starts  with  some  set  of  points  whose 
relative  positione  are  known  and  gradually  attaches 
new  sets  of  points  whoso  relative  poaitions  are 
computed  with  respeot  to  the  original  nucleus. 
Such  a process  may  be  viewed  as  a (possibly 
parallel)  solidification  of  parts  of  the 
measurement  graph  into  bodies. 

To  be  able  to  describe  incremental  construction 
processes  we  need  to  introduoe  a suitable 
formalism.  In  the  following  we  shall  describe  such 
a formalism,  then  apply  it  to  develop  anti  implement 
construct  ion  algorithm*. 

3.2.1  Stiff  hypergraphs 

A iivoergyaph  H ~ <V,2>  consists  of  a set  of 
vertioea  V and  a set  of  edges  2.  An  edge  is  a 
subset  of  vertioes  (not  necessarily  Just  two,  as  in 
graphs).  We  shall  use  this  generalised  notion  of 
an  edge  to  describe  a set  of  vertioes  whose 
relative  poaitions  are  known.  An  edge  is  said  to 
be  incident  upon  a given  vertex  if  it  contains  the 
vertex.  A vertex  is  said  to  be  iaoident  upon  a 
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figure  12.  Constructing  flexible  graphs  witn 
arbitrary  node  degrees 


figure  13.  A o-cooneoted  flexible  graph 


figure  lb,  A 5-cooaaated  flexible  graph 


given  edge  If  it  is  contained  in  this  edge.  We 
shall  consider  only  hypergraphs  for  which  the 
intersection  of  any  two  sd;es  contains  at  most  one 
point. 

To  develop  a visual  intuition  of  the  prooesses  to 
be  discussed,  one  should  insider  hypergraphs  as  a 
generalization  of  structures.  Rather  than 
considering  pin- Jointed  Pars,  we  consider 

pin- jointed  metal  sheets  of  an  arbitrary  shape  and 
nuaber  of  pins  (each  metal  sheet  corresponding  to 
an  edge  of  the  hypergraph  and  each  vertex 

corresponding  to  a Joint).  To  draw  further  on  the 
analogy,  we  shall  use  the  tern  link  as  an 

alternative  to  an  edge;  the  term  joint  will  be 
employed  as  an  alternative  to  vertex. 

Lot  d,  denote  the  degree  of  the  1-th  vertex  (l.e., 
the  nuaber  of  edges  incident  upon  it ) . Let  d^ 

denote  the  degree  of  the  i-th  edge  (l.e.,  the 
number  of  vertices  incident  upon  it).  The  degree 
of  tne  hypergrapn  H is  defined  to  be  the  number 

i(H)  * d,  , wnere  a(n)  is  the  total  oua- 

!»/  (./  * 

her  of  edges  (vert toes). 

The  nunoer  f(N)  = an^o-adCH)  is  called  the 
internal,  fcsadaa  ai  a,.  it  is  easy  to  verify  that 
for  a graph  H the  definition  of  f(H)  above 
degenerates  to  the  nuaber  of  internal  degrees  of 
freedom  defined  la  previous  aeotioaa. 

The  notion  of  stiffness  asn  be  generalised  to 
nypergcspas  as  follows,  a hypergreph  H is  said  to 
oe  stiff  iff  it  contains  a spanning  hypergreph  H* 
sunn  that 

1.  f(H' ) « 0 

2.  for  any  subhypergraph  U*  of  H'  f(U*)  > 0 

Let  js  introduce  a partial  order  over  hypergnpha, 
union  we  call  jaiiliait*  * 'Ypergreph  H*  Is  said  to 
be  a welding  of  a nypergrapfl  if  each  edge  of  is 
a union  of  edges  in  U which  span  a stiff 

lUbhypergrapti  of  h. 

It  is  possible  to  snow  that  stiffness  la  preserved 
under  welding.  Moreover,  each  nypertrapn  possesses 
e caique  weldiag  which  Is  aaxlnal  (cannot  be  welded 
any  core).  If  M(H)  designates  the  maximal  welding 
of  the  hypergreph  J),  then  f(M(K))  defir.ea  the 
deareea  it  frsedca  of  H.  Tt  can  be  shewn  that 
f<MtH ) ) > 3 with  equality  iff  H Is  stiff,  la  which 
case  H(H)  has  » single  edge  covering  all  the  nodes. 

3.2.2  icoresental  construct  ion  algorithms 

An  incremental  construction  algorithm  s a process 
that  starts  with  a given  positioning  problem  and 
develops  a solution  by  gradually  increasing  the 
sets  of  points  wnote  relative  locations  are  known. 
At  each  stage,  the  state  of  the  oasputatlon  may  be 
described  as  a hypergraph  whose  edges  consist  of 
sets  of  poihts  wnoee  relative  positions  ere  already 
known.  Such  a sypergrapn  Is  necessarily  a welding 


of  hypergraphs  in  the  previous  stages.  In  short, 
an  incremental  construction  algorithm  is  a process 
that  traces  a chain  in  the  partial  order  of 
welding. 

We  define  a welder  to  be  an  operator  that  takes  a 
hypergraph  and  produces  a 'welding  of  it . An 
Incremental  construction  algorithm  is  thus  a 
process  of  successive  applications  of  welders. 

We  have  developed  software  to  represent  and 
manipulate  structures  and  hypergraphs.  Two  types 
of  welders  have  been  implemented  and  some  simple 
construction  algorithms  tried.  We  possess  the 
tools  which  are  necessary  to  davalop 
position-locating  algorithms  of  increasing 
sophistication. 
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